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Abstract 

Due  to  the  unusually  high  heats  of  vaporization  of  room-temperature  ionic  liquids  (RTILs), 
volatilization  of  RTILs  through  thermal  decomposition  and  vaporization  of  the  decomposition 
products  can  be  significant.  Upon  heating  of  these  RTILs  in  vacuum,  their  gaseous  products 
were  detected  experimentally  via  tunable  vacuum  ultraviolet  photoionization  mass  spectrometry 
performed  at  the  Chemical  Dynamics  Beamline  9.0.2  at  the  Advanced  Light  Source. 
Experimental  evidence  for  di-  and  tri-alkylimidazolium  cations  and  nitrile  containing  anions 
indicates  thermal  decomposition  occurs  through  two  pathways:  deprotonation  of  the  cation  by 
the  anion,  and  dealkylation  of  the  imidazolium  cation  by  the  anion.  Additional  evidence 
supporting  these  mechanisms  was  obtained  using  thermal  gravimetric  analysis/mass 
spectrometry,  gas  chromatography/mass  spectrometry  and  temperature-jump  infrared 
spectroscopy.  In  order  to  predict  the  overall  thermal  stability  in  these  ionic  liquids,  the  ability  to 
accurately  calculate  both  the  basicity  of  the  anions  and  their  nucleophilicity  in  the  ionic  liquid  is 
critical.  Both  gas  phase  and  condensed  phase  (polarizable  continuum  model)  density  functional 
theory  calculations  provide  support  for  the  decomposition  mechanisms  and  a  means  to  determine 
thermal  stabilities  for  ionic  liquids  in  general. 

Keywords:  ionic  liquids,  thermal  decomposition,  synchrotron,  VUV,  photoionization  efficiency, 
polarizable  continuum  model. 
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INTRODUCTION 


Ionic  liquids  (with  m.p.  <100  °C)  consist  of  discrete  ions  and  have  distinctly  different 
physieal  and  chemieal  properties  from  moleeular  liquids.  The  most  widely  reeognized 
eharaeteristie  of  ionie  liquids  is  their  negligible  vapor  pressure,  making  them  attractive 
eandidates  for  replacements  of  volatile  organic  compounds  commonly  used  as  reaetion  solvents. 
In  order  to  depress  the  melting  point  of  ionic  solids,  it  is  necessary  to  reduee  the  Coulombie 
attraetion  between  cations  and  anions.  This  can  be  achieved  through  funetionalization  of  small, 
tightly-bound  ions  to  ereate  larger,  more  diffusely  charged  ions.  Often,  the  reduction  in 
Coulombie  attraction  between  the  ions  results  in  glass  formation  rather  than  aetual  melting  of  the 
solid.  In  the  suecessful  synthesis  of  ionic  liquids,  their  viseosities  tend  to  remain  unusually  high, 
which  can  limit  their  usefulness  as  solvents. 

The  dicyanamide  anion-based  family  of  room-temperature  ionic  liquids  (RTILs)  exhibits 
some  of  the  lowest  viscosities  of  any  known  RTILs.'  They  can  be  easily  prepared  from  the 
appropriate  halide  salts  and  silver  dieyanamide  and  were  the  first  RTILs  to  demonstrate 
hypergolie  reactivity  (spontaneous  ignition)  when  mixed  with  white  fuming  nitrie  aeid.  ’  Ionic 
liquids  are  not  only  well  known  for  their  negligible  vapor  pressures  but  also  for  their  thermal 
stability.  Because  of  their  extremely  low  affinity  for  protons,  bis-trifluoromethylsulfonylamide 
(NTf2‘)  based  RTILs  can  be  thermally  stable  up  to  almost  700  K.  While  significant  theoretical 
work  has  been  done  on  RTILs  in  the  gas  phase  and  by  force  field  methods  and  quantum 

o  n 

mechanics/moleeular  meehanies  (QM/MM),  '  predicting  reaction  energetics  in  the  condensed 

phase  is  diffieult  due  to  the  ionie  nature  of  the  solvent;  the  solvent  effeets  on  the  relative 

stabilities  of  the  reaetants,  transition  states  and  products  involved  in  the  reaction  will  effect  both 

the  rate  and  free  energy  of  reaetion.  The  development  and  implementation  of  an  aceurate  yet 
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efficient  theoretical  approach  to  investigate  the  condensed  phase  reaction  energetics  of  ionic 
liquids  is  highly  desirable.  Recent  work  in  this  regard  using  the  polarizable  continuum  model 

O 

(PCM)  shows  great  promise  in  this  application. 

In  this  paper,  we  investigate  the  thermal  stability  of  RTILs  with  CN  containing  anions 
such  as  dicyanamide  using  a  combination  of  several  experimental  probes  and  theory.  Thermal 
gravimetric  analysis  (TGA)  yields  the  enthalpies  of  vaporization,  and  when  coupled  with 
electron-impact  ionization  mass  spectrometry,  the  activation  enthalpies  for  individual  thermal 
decomposition  pathways  of  the  ionic  liquid  can  be  obtained.  Identification  of  vaporized  products 
is  also  confirmed  using  gas  chromatography/mass  spectrometry  (GC/MS),  thermal  gravimetric 
analysis-Fourier  transform  infrared  (TGA-FTIR)  spectroscopy  and  temperature-jump  (T-jump) 
FTIR  techniques.  Finally,  tunable  vacuum  ultraviolet  photoionization  time  of  flight  mass 
spectrometry  (VUV-PI-TOFMS)  is  used  to  identify  the  masses  and  ionization  potentials  of  the 
gas-phase  species  evolved  above  a  heated  RTIL  in  a  vacuum.  The  trends  in  these  results  have 
led  to  predictive  capabilities  using  density  functional  theory  (DFT)  and  direct  dynamics 
modeling  (see  accompanying  paper)  that  should  be  applicable  to  predicting  thermal  stabilities  of 
other  RTILs. 

EXPERIMENTAL 

I-ButyI-3-methyIimidazlium  tricyanomethanide  (BMIM^TCM',  >98%  purity)  and  I- 

butyI-3-methyIimidazolium  dicyanamide  (BMIM^dca',  >99%  purity)  were  purchased  from 

Merck  and  used  without  further  purification.  I  -Ethyl-3 -methylimidazolium  dicyanamide 

(EMIM'^dca'  >98.5%  purity),  1-ethylimidazolium  thiocyanate  (EMIM^SCN',  >99%  purity)  and 

1 -butyl-3 -methylimidazolium  thiocyanate  (BMIM^SCN',  >95%  purity)  were  purchased  from 

Fluka  and  used  without  further  purification.  1 -Butyl-3 -methylimidazolium  vinylogous 
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dicyanamide  (BMIM^vdca  )  was  synthesized  as  described  in  reference  l-Ethyl-2,3- 
dimethylimidazolium  dicyanamide  (EMMIM'^dca')  was  synthesized  according  to  the  procedure 
in  reference 

The  T-jump  experiment*'^’"  was  carried  out  on  a  rapid  scan  FTIR  spectrometer  with  4  cm' 
*  resolution  and  10  msec  time  resolution  between  spectra.  An  ionic  liquid  sample  of  about  3  pL 
was  analyzed  after  a  IF  capacitor  discharged  approximately  17  volts  across  a  nichrome  ribbon, 
reaching  maximum  temperature  in  ~50  msec.  The  discharge  of  the  capacitor  triggered  the  FTIR 
spectrometer.  The  experiment  was  conducted  with  an  initial  temperature  of  298  K  and  at  1 
atmosphere  pressure.  Maximum  temperatures  recorded  for  the  IFs  are  as  follows:  BMIM^SCN' 
Tmax  =  461  K,  BMlM+dca'  T^ax  =  684  K  and  BMIM+TCM'  T^ax  =  675  K. 

12  13 

The  TGA  and  TGA-MS  experiments  were  similar  to  those  described  previously.  ’ 
Briefly,  the  non-isothermal  TGA  measurements  were  carried  out  with  a  temperature  ramp  rate  of 
10  K/min  from  303-973  K,  and  the  isothermal  TGA  measurements  typically  were  carried  out 
from  423-573  K  in  either  5  K  or  10  K  steps  with  a  30  minute  hold  time  at  each  step. 

Gas  chromatography-mass  spectrometry  (GC-MS)  analysis  was  run  with  the  electron 
impact  ionizer  tuned  to  15  eV.  VUV-Pl-TOFMS  experiments  included  the  improved  effusive 
ionic  liquid  source*^  and  the  aerosol  source*"*  described  previously. 

THEORY 

Hybrid  density  functional  theory  was  applied  at  the  M06/6-31+G(d,p)  level  of  theory*^  to 

investigate  these  IF  systems.  M05  has  been  shown  to  provide  a  good  level  of  theory  for 

calculating  energetics  of  ionic  liquid  systems  (mean  average  deviation,  MAD,  of  9.5  kJ/mol),*^ 

and  M06  is  presumed  to  be  as  good  if  not  better  than  M05  for  our  purposes.  Stationary  states  of 

reactants,  intermediates,  transition  states  and  products,  as  well  as  adiabatic  photoionization 
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potentials  were  calculated  using  Gaussian09.  Transition  states  were  confirmed  to  have  a  single 
imaginary  frequency  and  through  internal  reaction  coordinate  analysis. 

A  variant  of  the  quantum  mechanical  continuum  solvation  model  SMD  known  as  the 
generic  ionic  liquid  (SMD-GIL)  model'^  was  employed  to  calculate  free  energies  in  the  ionic 
liquid  condensed  phase.  The  SMD-GIL  solvent  descriptor  input  parameters  used  for  EMIM^dca' 
as  a  solvent  are  e  =  11.50,  n  =  1.4300,  y  =  61.24,  =  0.229, 1.^2  ^  0.265,  <p=  0.2308,  and  \j/= 

0.0000,'*^  while  (p  =  0.2727  and  0.200  for  EMIM^SCN'  and  EMIM^TCM',  respectively.  The 
fraction  of  non-hydrogen  atoms  that  are  electronegative  halogen  atoms  ( is  zero  for  dca',  SCN' 

,  TCM’  and  vdca’  containing  ionic  liquids  and  y/=  0.0500  for  EMIM^Bf.  SMD-GIL 

1  7 

calculations  were  carried  out  using  M06/6-31+G(d,p)  in  Gaussian  09. 

RESULTS 

The  thermal  decomposition  of  CN  containing  anionic  RTILs  was  investigated  using 
several  complementary  experimental  probes,  including  thermal  gravimetric  analysis,  thermal 
gravimetric  analysis  coupled  to  either  a  mass  spectrometer  or  a  Fourier  transform  infrared 
(FTIR)  spectrometer,  temperature-jump  FTIR  spectroscopy,  GC/MS  and  tunable  vacuum 
ultraviolet  photoionization  time  of  flight  mass  spectrometry.  For  brevity,  a  representative 
spectrum  from  each  experimental  technique  will  be  presented  here  and  the  remainder  of  the  data 
can  be  found  in  the  supporting  information.  Numeric  data  is  summarized  in  Tables  1-5. 

Thermal  Gravimetric  Analysis 

Thermal  gravimetric  analysis  (TGA)  is  a  measurement  of  the  mass  of  a  sample  as  it  is 
heated,  typically  with  an  inert  gas  flowing  over  the  sample.  In  an  inert  environment,  the  mass 
loss  as  the  temperature  increases  can  be  the  result  of  vaporization  of  the  sample  or  of  the  thermal 


decomposition  of  the  sample  followed  by  subsequent  vaporization  of  the  thermal  decomposition 
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products.  For  volatile  species,  simple  evaporation  takes  place  and  100%  of  the  sample  mass  is 
lost  due  to  vaporization  of  the  intact  species.  The  rate  of  evaporation  is  dependent  on  the 
enthalpy  of  vaporization  of  the  species,  the  surface  area  of  the  sample,  and  the  flow  and 
molecular  mass  of  inert  gas  over  the  sample.  ’  For  easily  decomposing  species,  first 
decomposition  must  occur  and  the  products  must  diffuse  to  the  surface  of  the  sample,  where  they 
can  evaporate.  Mass  loss  due  to  thermal  decomposition  processes  can  be  less  than  100%  if 
thermal  decomposition  products  have  low  volatility  or  form  solids.  The  rate  of  thermal 
decomposition  depends  on  the  activation  enthalpies  for  all  decomposition  pathways  and  is  not 
dependent  upon  surface  area  of  the  sample.  ’  The  total  mass  loss  rate  by  TGA  is  the  sum  of 
the  mass  loss  from  evaporation  and  the  mass  losses  from  all  thermal  decomposition  processes 
which  produce  volatile  species.  However,  because  the  evaporation  and  thermal  decomposition 
rates  are  dependent  on  different  factors,  the  enthalpy  of  vaporization  for  low-volatility  species 
such  as  RTILs  can  be  accurately  determined  by  an  isothermal  analysis  and  the  thermal 
decomposition  rates  are  better  measured  by  fast  heating  rates.  When  the  enthalpy  of 
vaporization  is  comparable  to  or  higher  than  the  lowest  activation  enthalpy  for  thermal 
decomposition,  then  mass  loss  by  thermal  decomposition  competes  with  or  dominates  over  mass 
loss  by  evaporation.  ’  Enthalpies  of  vaporization  of  RTILs  range  from  120  kJ/mol  to  over  200 
kJ/mol,  and  thermal  decomposition  is  a  highly  competitive  process  in  these  compounds. 

The  TGA  of  l-ethyl-3-methylimidazolium  dicyanamide,  EMIM^dca',  indicates  a  mass 

loss  onset  (non-zero  slope)  temperature  above  510  K  with  two  distinct  peaks  at  590  and  690  K 

(Figure  la).  Fitting  the  lowest  temperature  peak  in  an  Arrhenius  plot  (Figure  lb)  yields  an 

enthalpy  for  the  overall  mass  loss  rate  of  151.7  ±  14.5  kJ/mol  (±  95%  confidence  limits)  at  an 

average  temperature  (Tavg)  of  566  K  (Table  1).  In  order  to  determine  the  enthalpy  of 
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_i_  1  o  22 

vaporization  (A/fyap)  of  EMIM  dca',  isothermal  TGA  analysis  (described  previously)  ’  yields  a 
A//vap  value  of  A//vap  =  128.3  ±  19.2  kJ/mol  (Tavg  =  473  K),  which  when  corrected  to  298  K  using 
Ag  Cp  of  -72.0  J/K  (see  Table  1  and  discussion  for  details)/  gives  Af/vap(298)  =  140.9  ±  19.2 
kJ/mol. 


1/T 


Figure  1.  Thermal  gravimetric  analysis  of  1 -ethyl-3 -methylimidazolium  dicyanamide:  a)  TGA, 
b)  Arrhenius  plot  of  slope  from  Figure  la,  c)  isothermal  TGA  to  determine  Af/yap  =  140.9  ±  19.2 
kJ/mol  (298  K). 

Table  1.  Experimental  enthalpies  of  activation,  A/7^ (overall),  and  vaporization,  A/fyap  (at  Tavg 
and  298  K)  determined  by  TGA.  Enthalpies  are  in  kJ/mol,  Af  is  in  J/mol»K  and  temperatures 
are  in  K. 


Ionic  liquid 

AEf*(overall) 

Tavg 

AHyap 

Tavg 

ASfO 
^1  '^pm 

A^vaD(298) 

EMIM^dca 

151.7  ±  14.5 

566 

128.3  ±  19.2 

473 

-72.0“  ± 

3.3 

140.9  ±  19.2 

BMIM^dca 

159.5  ±  10.2 

560 

133.1  ±5.0 

488 

-72.5“  ± 

3.3 

146.9  ±5.0 

EMMIM^dca 

268.4  ±50.9 

618 

185.2  ±  14.2 

513 

-70 

200.3  ±  14.2 

EMIM^SCN' 

118.0  ±3.7 

535 

150.7  ±  10.8 

458 

-H 

J=> 

00 

o 

2.8 

160.4  ±  10.8 

BMIM^SCN' 

165.7  ±5.0 

517 

168.6  ±9.2 

448 

-H 

j=> 

o 

3.2 

179.9  ±9.2 

BMIM^TCM’*^ 

183.0  ±4.3 

582 

126.0  ±6.2 

473 

-78.9"  ± 

3.6 

139.8  ±6.2 

BMIM^dca 

105.1  ±25.3 

475 

62.1  ±25.3 

453 

-70 

73.0  ±25.3 

“reference 

^reference  ^references  ^estimated  value-  see  text  i 

'or  details. 

For  EMlM^dca'  and  several  other  CN  containing  anionic  RTlLs,  there  are  several  peaks 
in  the  TGA  and  the  mass  loss  does  not  reach  100%  at  high  temperature.  This  indicates  that  there 
are  several  different  thermal  decomposition  processes  occurring  in  different  temperature 
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regimes,  and  that  some  solid  is  formed  which  is  not  volatile.  Previous  work  on  dicyanamide- 
based  ionic  liquids  indicates  that  the  dicyanamide  anion  can  polymerize  when  heated,  perhaps 

9Q  ^0 

through  a  melamine  intermediate.  ’  This  will  be  discussed  in  more  detail  in  the  discussion 
section. 

Thermal  Gravimetric  Analysis-Mass  Spectrometry 

A  powerful  enhancement  to  standard  TGA  is  to  couple  the  TGA  to  a  mass  spectrometer 
via  a  heated  capillary  tube.  Volatile  species  evolved  from  the  TGA  can  then  be  identified  by 
their  mass-to-charge  ratios  (m/z)  and  their  fragmentation  patterns.  Furthermore,  by  monitoring  a 
product  species  parent  mass  as  a  function  of  temperature,  the  activation  enthalpy  for  its 
formation  can  be  calculated  using  an  Arrhenius-type  analysis.  The  major  species  detected  by 
TGA-mass  spectrometric  analysis  of  EMIM^dca'  are  m/z  =  82  and  96  (Figure  2a)  and  an 
Arrhenius  plot  for  these  masses  (Figure  2b)  yields  the  activation  enthalpies  for  their  formation, 
A//^96  =  160  ±  16  kJ/mol  and  175  i  17  kj/ixiol  3.1  Tg^yg  505  Iv.  .A.  coiiS0rv3tiv0  0Stim3t0 

of  ±  10%  was  assigned  to  the  uncertainties  in  A//^  due  to  possible  systematic  errors  arising  from 
mass  transport  (related  to  the  chemical’s  vapor  pressure)  differences  between  the  TGA  and  the 
mass  spectrometer  and  ion  transmission  efficiency  differences  in  the  mass  spectrometer  itself 
A/f^  values  calculated  for  selected  masses  of  all  ionic  liquids  are  seen  in  Table  2. 
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Figure  2.  a)  TGA-MS  analysis  of  EMIM^dca':  mass  82  (methylimidazole)  and  96 
(ethylimidazole),  and  b)  Arrhenius  plot  of  data  from  Figure  2a:  A//^96  =  159.6  ±  2.3  kJ/mol  (505 
K),  Aff^82  =  174.8  ±  2.4  kJ/mol  (505  K). 


Table  2.  Experimental  enthalpies  of  activation,  A//^,  for  individual  masses  (at  Tavg)  determined 
by  TGA-MS  analysis.  Enthalpies  are  in  kJ/mol,  m/z  values  are  in  amu  and  temperatures  are  in 
K.  Uncertainties  are  estimated  to  be  ±  10%. 


Ionic  liquid 

m/z 

m/z 

T 

Aavg 

EMIM+dca 

160  ±  16 

96 

175  ±  17 

82 

505 

BMIM+dca 

150  ±  15 

82 

162  ±  16 

97 

510 

EMMIM+dca 

140  ±  14 

96 

137  ±  14 

no 

530 

BMIM+TCM' 

128  ±  13 

78 

98±  10 

82 

525 

EMIM+SCN' 

128  ±  13 

73 

109  ±  11 

82 

478 

126  ±  13 

87 

110±  11 

96 

478 

Vacuum  Ultraviolet  Photoionization  Time  of  Flight  Mass  Spectrometry 

Representative  vacuum  ultraviolet  time-of- flight  mass  spectrometry  (VUV-TOFMS)  data 
for  EMIM'^dca'  can  be  seen  in  Figure  3.  Figure  3a  shows  a  contour  plot  of  the  photoion 
intensities  with  a  VUV  photon  energy  range  of  7.4-10.0  eV  at  an  effusive  source  temperature  of 
T=473  K,  and  Figure  3b  is  the  corresponding  photoionization  mass  spectrum  at  9.0  eV. 
Photoionization  efficiency  (PIE)  curves  for  masses  82,  96,  110  and  150  are  seen  in  Figure  3c-d, 
and  the  total  ion  intensities  (sum  over  all  photon  energies)  of  these  masses  are  plotted  as  a 
function  of  effusive  source  temperature  in  Figure  3e.  Photoion  appearance  energies  were 
determined  for  the  major  photoions  by  fitting  the  PIE  curves  to  the  function  a(E-Eo)^  near  the 
photoionization  threshold,  as  described  previously,"^  '"^’^'  and  are  reported  along  with  their 
calculated  (M06/6-31+G(d,p),  298  K)  adiabatic  ionization  energies  and  relative  ion  abundances 
in  Table  3. 
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a 


Figure  3.  a)  VUV-PI-TOFMS  analysis  of  EMIM+dca'  at  473  K  from  7.4  to  10.0  eV  photon 
energy,  and  b)  9.0  eV  photon  energy,  c)  PIE  curves  for  masses  in  (b),  d)  expanded  view  of  the 
photoion  thresholds  from  (c),  and  e),  total  ion  count/relative  abundance  of  the  photoions  from 


450  to  525  K. 
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Table  3.  VUV-TOFMS  data  for  all  ionic  liquids  investigated.  AG298  values  are  in  eV. 


EMmrdca'.  463  K 

mass  (amu) 

appearance  energj'  (e\^ 

AG2geM06/6-31±G(d,p) 

relative  abundance  (%) 

27 

13.4  ±0.2 

13.3  ±0.2 

4.4 

82 

8.3  ±  0.2 

8.4  ±0.2 

14.2 

96 

8.3  ±  0.2 

8.3  ±0.2 

23.9 

110 

7.6  ±  0.2 

7.7  ±0.2 

11.4 

150 

7.5  ±0.2 

7.6  ±0.2 

45.3 

BMnTdca'.  483  K 

mass  (amu) 

appearance  energj'  (e\^ 

AG29eM06/6-31±G(d,p) 

relative  abundance  (%) 

27 

13.4  ±0.2 

13.3  ±0.2 

4.0 

82 

8.4  ±  0.2 

8.4  ±0.2 

20.5 

97 

7.6  ±  0.2 

ND 

18.0 

124 

8.2  ±  0.2 

8.3  ±0.2 

15.6 

137 

7.9  ±  0.2 

ND 

15.1 

178 

7.6  ±  0.2 

8.3  ±0.2 

5.4 

EMMnTdca'.  423  K 

mass  (amu) 

appearance  energj'  (e\^ 

AG29eM06/6-31±G(d,p) 

relative  abundance  (%) 

96 

7.7  ±  0.2 

8.0  ±0.2 

17.0 

110 

7.9  ±  0.2 

7.9  ±0.2 

6.1 

124 

6.0  ±  0.5 

6.1  ±0.2 

52.8 

125 

7.2  ±  0.5 

7.2  ±0.2 

23.9 

EMIM'SCN',  467  K 

mass  (amu) 

appearance  energj'  (e\^ 

AG  29g  M06/6-31±G(d,p) 

relative  abundance  (%) 

27 

13.3  ±0.2 

13.3  ±0.2 

3.4 

73 

8.2  ±0.2 

9.1  ±0.2(CH3NCS) 

23.2 

82 

8.4  ±  0.2 

8.4  ±0.2 

8.8 

87 

9.1  ±0.2 

9.0  ±0.2(CH3CH2NCS) 

4.6 

96 

8.3  ±  0.2 

8.3  ±0.2 

30.9 

142 

7.1  ±0.3 

7.2  ±0.2 

29.1 

BMarSCN",  493  K 

mass  (amu) 

appearance  energj'  (e\^ 

AG29gM06/6-31±G(d,p) 

relative  abundance  (%) 

27 

13.4  ±0.2 

13.3  ±0.2 

4.4 

73 

9.1  ±0.2 

9.1  ±0.2(CH3NCS) 

29.3 

82 

8.3  ±  0.2 

8.4  ±0.2 

19.7 

97 

9.3  ±  0.2 

ND 

9.3 

124 

8.1  ±0.2 

8.3  ±0.2 

15.1 

137 

9.2  ±  0.2 

ND 

20.3 

170 

8.8  ±0.2 

7.1  ±0.2 

1.8 

BMnrrcM'.  433  k 

mass  (amu) 

appearance  energj'  (e\^ 

AG2sgM06/6-31±G(d,p) 

relative  abundance  (%) 

82 

8.2  ±  0.2 

8.4  ±0.2 

1.1 

96 

7.6  ±  0.2 

ND 

67.9 

110 

7.7  ±  0.2 

ND 

2.0 

124 

8.2  ±0.2 

8.3  ±0.2 

4.6 

137 

7.6  ±  0.2 

ND 

17.8 

139 

7.1  ±0.3 

6.9  ±0.2 

3.4 

229 

6.6  ±  0.5 

7.3  ±0.2 

3.1 

BMnrvdca\413K 

mass  (amu) 

appearance  energj'  (e\^ 

AG  29g  M06/6-31±G(d,p) 

relative  abundance  (%) 

96 

7.7  ±  0.2 

ND 

1.2 

110 

8.0  ±  0.2 

ND 

8.0 

124 

8.2  ±  0.2 

8.3  ±0.2 

2.1 

137 

7.6  ±  0.2 

ND 

8.0 

178 

7.1  ±0.3 

8.3  ±0.2 

43.6 
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VUV-PI-TOFMS  for  1 -butyl-3 -methylimidazolium  dicyanamide  (BMIM'^dca  )  was 
extensively  performed  for  this  species.  Figures  4a-c  represent  two  VUV-TOFMS  spectra  which 
are  generated  by  low  temperature  evaporation  of  BMIM^dca'  aerosols  as  described  previously^"^ 
(Figures  4a  and  b,  380  K),  or  by  the  effusive  ionic  liquid  source  at  483  K  and  photoionized  at 
10.0  eV  (Figure  4c).  Figure  4d  is  from  gas  chromatograph-mass  spectrometric  (GC/MS) 
analysis  of  BMIM^dca'  with  a  15  eV  electron-impact  ionization  source,  indicating  a  parent  ion 
mass  of  m/z  =178.  In  Figure  5,  the  scaled  mass  spectra  indicate  the  changes  in  the  BMIM"^ 
cation  signal  as  a  function  of  temperature  in  the  range  of  438-483  K  photoionized  at  12.0  eV  and 
averaged  over  5  million  repeller  pulses. 


Figure  4.  a)  VUV-PI-TOFMS  spectrum  of  BMIM^dca'  aerosol  at  380  K  and  10.0  eV  and  b) 
expanded  view  of  the  BMIM^  cation  at  m/z  139  from  (a),  c)  BMIM^dca’  spectrum  from  the 
effusive  source  at  483  K  and  10.0  eV  and  d)  GC/MS  analysis  of  BMIM'^dca'  with  15  eV  electron 
impact  ionization  energy. 
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Figure  5.  VUV-PI-TOFMS  spectram  of  BMIM^dca'  for  the  effusive  souree  in  the  range  of  438- 
483  K  at  12.0  eV  photoionization  energy,  5  million  shots,  traces  scaled  for  comparison. 

Temperature-jump  Fourier  Transform  Infrared  Spectroscopy 

Likely  due  to  polymerization,  no  coproducts  containing  the  dca’  or  TCM'  anions  were 
observed  using  the  TGA-MS  or  the  effusive  source  in  the  VLFV  experiments.  Therefore, 
temperature-jump  Fourier  transform  infrared  spectroscopy  (T-jump  FTIR)  was  employed  in  an 
attempt  to  rapidly  vaporize  the  products  before  they  could  polymerize.  The  experimental  design 
is  based  on  previous  work  by  Brill  and  coworkers. T-jump  FTIR  spectroscopy  results  for 
BMIM^SCN',  BMIM^TCM’  and  BMIM^dca’  are  seen  in  Figure  6. 
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Figure  6.  T-jump  FTIR  spectra  of  a)  BMIM^SCN'  decomposition  products  (black)  compared  to 
a  H3CNCS  (blue,  ref  and  3:1  ratio  of  butylimidazole  and  methylimidazole  (red),  b) 
BMIlS/f^TCM'  decomposition  products  (black)  and  linear  ketenimine  (red,  from  reference  ^^),  c) 
an  expanded  view  of  the  central  peak  in  (b),  d)  BMIM'^dca’  decomposition  products  (black)  and 
butylimidazole  (red),  e)  DFT  IR  simulation  of  methylated  TCM',  and  f)  DFT  IR  simulation  of 
methylated  dca'. 

DISCUSSION 

Thermal  Gravimetric  Analysis 

In  comparing  the  TGA  data  for  all  of  the  investigated  CN-containing  anionic  ionic 
liquids,  the  only  IL  to  achieve  near  100%  mass  loss  are  the  thiocyanates  l-ethyl-3- 
methylimidazilium  thiocyanate,  EMIM^SCN',  and  1 -butyl-3 -methylimidazolium  thiocyanate, 
BMIlVf^SCN'.  All  of  the  other  ILs  investigated  have  mass  losses  <100%,  indicating  solid  residue 
formation,  likely  through  polymerization  of  the  anions  or  their  thermal  decomposition 
products.^^’^^  This  indicates  that  the  thermal  decomposition  products  in  the  SCN’  systems  either 
do  not  easily  polymerize  or  are  able  to  vaporize  before  polymerization  can  occur.  From  Table  1, 
based  on  the  overall  values,  the  trend  in  stability,  from  least  to  most  stable,  is:  BMIM^vdca’ 
DISTRIBUTION  A:  Approved  for  public  release,  distribution  unlimited 


14 


<  EMIM+SCN'  <  EMIM+dca  <  BMIM+dca  <  BMIM+SCN'  <  BMIM+TCM'  <  EMMIM+dca. 


Similarly,  from  Table  1,  based  on  the  overall  Aifyap  values,  the  trend  in  enthalpies  of 
vaporization,  from  lowest  to  highest,  is:  BMIM^vdca'  <  BMIM^TCM’  <  EMIM^dea'  < 
BMIM^dea’  <  EMIM^SCN’  <  BMIM^SCN'  <  EMMIM^dca’,  whieh  is  consistent  with  previous 
experimental  and  theoretical  heats  of  vaporization  for  these  The  heat  of  vaporization 

determined  for  BMINT^vdca'  is  unusually  low  (Table  1),  and  is  likely  attributable  to  mass  loss 
from  thermal  decomposition  instead  of  from  vaporization  of  the  intact  IE.  If  this  is  true,  then 
based  on  enthalpies  of  vaporization,  the  most  volatile  IE  is  BMIM^TCM'.  This  species  has  been 
detected  directly  in  the  vapor  phase  and  has  one  of  the  highest  vapor  pressures  of  known  ILs.'^’^^ 
Clearly,  for  the  thiocyanate  ILs,  since  their  overall  activation  enthalpies  determined  by  TGA  are 
significantly  lower  than  their  enthalpies  of  vaporization,  thermal  decomposition  is  strongly 
favored  over  vaporization  at  moderate  temperatures  (T<573  K)  in  these  ILs.  Conversely,  for  the 
tricyanomethanide  and  trialkylimidizolium  dicyanamide  ILs,  their  enthalpies  of  vaporization  are 
significantly  lower  than  their  overall  thermal  decomposition  activation  enthalpies,  and 
vaporization  is  favored  over  thermal  decomposition  at  moderate  temperatures  (T<573  K),  as  has 
been  observed  previously  for  dialkylimidazolium  carboxylate  ionic  liquids.^^  For  the 
dialkylimidazolium  dicyanamides,  where  their  overall  thermal  decomposition  activation 
enthalpies  and  enthalpies  of  vaporization  are  fairly  similar,  it  is  likely  that  thermal  decomposition 
and  vaporization  are  competing  below  T  =  573  K.  It  is  important  to  note  that,  according  to 
Trouton’s  Rule  (ASvap  ~  10. 5R),  entropic  contributions  are  typically  larger  for  vaporization  than 
for  thermal  decomposition  processes  in  the  condensed  phase  discussed  here  (AS  <  4R),  and  this 
effectively  lowers  the  free  energy  of  vaporization  more  than  for  thermal  decomposition  versus 

considering  only  the  enthalpies  in  these  two  competing  processes. 
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In  order  to  accurately  determine  the  enthalpy  of  vaporization  of  ionic  liquids  under 
standard  conditions  (298  K),  it  is  important  to  accurately  determine  AfCpm,  the  change  in  molar 
heat  capacity  from  liquid  to  the  gas  phase.^^’^^’"^*^''^^  A  recent  study  by  Verevkin^  and  coworkers 
have  demonstrated  that  Af  Cp^  can  be  estimated  by: 

AfC°„=  -2I?-(C°„-C°„),  (1) 

where:  (C°m '  (2) 

and  ttp  is  the  thermal  expansion  coefficient,  Kj  is  the  isothermal  compressibility,  and  Vm  is  the 
molar  volume  of  the  ionic  liquid.  Prior  to  this  work,  Af  Cp^,  values  were  simply  estimated  to  be 
100  J/mol*K  and  this  assumption  was  the  main  source  of  error  in  measuring  standard  enthalpies 
of  vaporization  of  ionic  liquids.  For  this  work,  literature  values  for  ap  and  Kt  for  five  of  the 
ionic  liquids  were  used  to  calculate  the  AfCp^  values  in  Table  1.  The  100  J/mol*K  value  for 
Af  Cpm  used  previously  was  based  on  much  work  with  the  bis-trifluorosulfonylamide  anion- 
based  ionic  liquids,  which  have  large  molar  heat  capacities  (500-1000  J/mol»K  at  298  K^).  The 
ionic  liquids  investigated  here  have  smaller  molar  heat  capacities  (-300  J/mol*K)  and  would  be 
expected  to  have  smaller  Af  Cp^  values  as  well.  In  fact,  the  AfCpm  values  in  Table  1  average 
about  70  J/mol*K,  and  this  value  was  used  as  an  estimate  for  Af  Cp^  for  EMMIM^dca'  and  for 
BMIM^dca',  where  no  ap  or  Kt  values  were  available.  This  more  accurate  value  for  Af  Cp^, 
decreases  the  standard  enthalpies  of  vaporization  for  these  ionic  liquids  by  an  average  of  5.1 
kJ/mol  relative  to  those  using  the  AfCp^,  =  -100  J/mol*K  value. 

Thermal  Gravimetric  Analysis-Mass  Spectrometry 
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Since  the  TGA  results  do  not  discriminate  between  vaporization  and  the  thermal 
decomposition  pathways  available,  TGA-MS  is  better  suited  for  analyzing  specific  reaction 
pathways  by  monitoring  volatile  thermal  decomposition  products  as  a  function  of  temperature; 
vaporization  of  ILs  is  not  likely  to  be  detected,  as  the  ILs  will  condense  in  the  capillary  coupling 
the  TGA  to  the  mass  spectrometer.  Major  mass  peaks  that  were  detected  are  listed  in  Table  2 
along  with  their  activation  enthalpies  of  formation,  which  were  calculated  using  an  Arrhenius 
type  analysis  on  the  low  temperature  regime  of  the  TGA-MS  plots  in  Figure  2a  and  in  the 
supplementary  information. 

TGA-MS  results  for  the  thermally  unstable  BMIM^vdca'  were  complex  and  difficult  to 
interpret,  likely  due  to  multiple  decomposition  pathways,  and  will  not  be  discussed  here.  The 
next  most  sensitive  IL  to  thermal  decomposition,  EMIM^SCN',  evolves  masses  at  m/z  73,  82,  87 
and  96  with  activation  enthalpies  of  (Tavg  =  478  K)  A/f^vs  =  128  ±  13  kJ/mol,  A77^82  =  109  ±11 
kJ/mol,  A//^87  =  126  ±  13  kJ/mol  and  Af7^96  =  110  ±  11  kJ/mol  (uncertainty  estimated  at  10%), 
respectively,  and  a  small  signal  at  m/z  142  which  was  too  small  to  determine  the  thermal 
decomposition  activation  enthalpy.  M/z  27  was  difficult  to  detect  due  to  interference  from  the 
large  m/z  28  N2  carrier  gas  signal.  Note  that  these  activation  enthalpies  are  quite  close  to  the 
overall  thermal  decomposition  activation  enthalpy  determined  by  TGA  of  A/7^  =  114.9  ±  3.6 
kJ/mol.  It  is  likely  that  these  products  are  CH3NCS  (73),  methylimidazole  (82),  CH3CH2NCS 
(87)  and  ethylimidazole  (96)  formed  through  an  Sn2  type  reaction  whereby  the  SCN'  anion 
attacks  the  alkyl  groups  on  the  imidazolium  cation.  Mass  142  will  be  discussed  in  the  next 
section. 

For  BMIM'^dca’  and  EMIM'^dca’,  the  major  product  masses  detected  were  82  and  97,  and 

82  and  96,  respectively.  While  the  82  and  96  masses  are  likely  methyl-  and  ethylimidazole 
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formed  via  an  Sn2  reaction  similar  to  the  thiocyanate  ILs,  the  structure  of  mass  97  isn’t  clear.  A 
very  small  signal  at  m/z  124  was  detected  for  BMlM^dca'  that  corresponds  to  the  formation  of 
butylimidazole,  but  the  signal  was  too  small  to  determine  a  thermal  decomposition  activation 
enthalpy  for  this  species.  The  relatively  small  signal  for  butylimidazole  is  perhaps  due  to  the 
lower  vapor  pressure  of  butylimidazole  versus  methyl-  and  ethylimidazole.  For  EMMlM^dca', 
the  primary  masses  detected  were  110  and  96,  corresponding  to  the  Sn2  formation  of  l-ethyl-2- 
methylimidazole  and  1,2-dimethylimidazole. 

Finally,  for  BMIM^TCM',  masses  were  detected  via  TGA-MS  at  78,  82,  96,  107  and  124. 
While  masses  at  m/z  82  and  124  are  likely  methylimidazole  and  butylimidazole,  the  source  of 
m/z  78,  96  and  107  aren’t  clear,  and  further  investigation  is  warranted. 

Vacuum  Ultraviolet  Photoionization-Time  of  Flight  Mass  Spectrometry 

It  is  accepted  that  for  many  aprotic  ionic  liquids  they  evaporate  as  intact  ion  pairs:  C^A' 

(/)  ->  C"  A'  (g)  and  that  upon  ionization  of  these  gaseous  ion  pairs,  dissociation  to  the  cation 
and  the  radical  A  is  facile:  C^A'  C^A  +  e'  ->  +  A  +  t  Due  to  the  very  high  detection 

sensitivity  afforded  at  the  Chemical  Dynamics  Beamline  at  the  Advanced  Light  Source,  the  mass 
signals  resulting  from  the  ion  pairs  and  those  from  thermal  decomposition  products  or  even  trace 
contaminants  or  previous  IL  experiments  needed  to  be  carefully  differentiated.  In  many  cases, 
the  thermal  decomposition  products  are  detected  at  lower  source  temperatures  than,  or 
simultaneously  with,  the  ionic  liquid  vapors. 

A  common  feature  amongst  the  dialkylimidazolium  ionic  liquids  investigated  in  this  work 

was  the  significant  loss  of  1  amu  from  the  EMIM^  cations,  loss  of  2  amu  from  the  BMIM^ 

cations,  and  loss  of  1  amu  from  the  EMMIM"^  cation  to  a  lesser  extent.  A  second  trend 

throughout  these  ILs  was  the  appearance  of  a  mass  peak  in  between  the  mass  of  the  cation  and 
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the  ion  pair,  namely  the  mass  of  the  ion  pair  -  27  amu.  Finally,  as  was  mentioned  previously, 

there  is  strong  evidence  of  Sn2  dealkylation  of  the  imidazolium  by  the  anions.  Previous  work  on 

ILs  with  more  basic  anions  has  determined  that  neutralization  can  occur  through  proton  transfer 

from  C2  on  the  imidazolium  to  the  anion,  leading  to  the  formation  of  a  neutral  carbene"^^  '^^  and 

an  acid  (protonated  anion)  which  can  then  vaporize  and  be  detected,  unless  the  protonated  anion 

can  undergo  polymerization  to  form  a  solid.  In  the  case  of  EMIM'^dca',  the  carbene  is  detected  at 

m/z  110  with  a  photoionization  energy  of  7.6  ±  0.2  eV  (Figure  3c)  and  matches  well  with  the 

M06  calculated  ionization  energy  for  the  EMIM  carbene  (EMIM:)  of  7.7  ±  0.2  eV.  While  m/z 

110  was  not  detected  in  the  case  of  EMUS/f^SCN',  the  formation  of  m/z  142  and  HCN  strongly 

suggest  that  these  products  are  formed  via  the  EMIM:  carbene  intermediate  (Figure  7).  The 

extent  of  carbene  formation  should  be  dependent  on  two  factors:  the  basicity  of  the  anion  and  the 

stability  of  the  carbene  itself  Table  4  illustrates  the  relative  strength  of  the  base  by  determining 

the  free  energy  of  acidity:  essentially  how  much  free  energy  it  takes  to  separate  a  proton  from  the 

acid  HA.  The  higher  the  value  of  AGacid,  the  more  stable  the  acid  HA  and  the  less  likely  the 

anion  will  remain  as  a  negative  anion  in  the  presence  of  acidic  hydrogens.  Species  with 

AGacid(g)  values  smaller  than  300  kcal/mol  (1255  kJ/mol)  are  considered  to  be  superacids"^^  and 

these  tend  to  make  good  candidates  as  anions  for  IE  synthesis.  The  generic  ionic  liquid  (SMD- 

GIL)  variant'^  of  the  polarizable  continuum  model  was  employed  here  to  calculate  the  AGacid(l) 

values  in  an  imidazolium  ionic  liquid  and  with  water  as  a  solvent  as  well  (Table  4).  The  solvent 

effect  on  these  systems  reduces  the  free  energy  AGacid  values  from  around  1200  kJ/mol  to  around 

520  kJ/mol  in  both  the  ionic  liquid  and  in  water.  This  reduction  in  energy  is  primarily  driven  by 

the  large  solvation  energy  of  a  proton  (~550  kJ/mol)  versus  anions  (~210  kJ/mol  on  average)  in 

the  condensed  phases  (table  in  SI).  Based  on  the  AGacid  values  in  Table  4,  the  expected  trend  for 
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formation  of  C^-1  amu  photoions  to  appear  in  the  VUV  mass  spectrum  would  be:  SCN'  >  dca'  > 
vdca'  >  TCM'.  With  the  exception  of  EMIM'^SCN'  which  has  a  short-lived  m/z  110  as  described 
below,  this  trend  is  confirmed  in  the  VUV  experiments,  and,  in  fact,  no  C^-1  amu  peak  was 
detected  for  BMIM^TCM'.  The  AGacid  values  for  loss  from  the  C2  of  EMIM^  and  BMIM^  or 
from  the  C2  methyl  of  EMMINT^  are  262.9,  264.6  and  261.8  kJ/mol,  respectively,  and  essentially 
are  the  same  within  the  margin  of  error  of  M06  and  therefore  the  anion  basicity  is  independent  of 
cation  for  these  cations.  Note  that  the  deprotonation  of  EMMIM^  leads  to  formation  of  a  stable 
methylene  product  at  m/z  124  and  not  a  reactive  carbene  species. 

The  nitrate  and  bromide  anions  were  included  in  this  analysis  to  compare  their  relative 
basicities  with  the  nitrile  containing  anions.  Previous  gas  phase  acidity  calculations  indicated 
that  the  dca’  anion  was  less  basic  than  the  NO3’  anion,  and  the  equilibrium  HNO3  +  dca’  NO3’ 
+  Hdca  would  favor  protonation  of  the  NOs’.^  From  these  SMD-GIL  calculations,  the  solvation 
effect  stabilizes  Hdca  more  than  HNO3  and  NO3’  more  than  dca’,  and  would  result  in  a 
spontaneous  proton  transfer  from  HNO3  to  dca’.  This  solvation  effect  might  help  to  explain  the 
high  reactivity  of  dicyanamide-based  ionic  liquids  towards  nitric  acid.^ 


Table  4.  Calculated  free  energies  of  acidity,  AGacid,  in  the  gas  phase  and  in  the  condensed  phase 


by  SMD-GIL  and  SMD  (water)  at  the  M06/6-31+G(d,p)  level  of  theory. 


AGacid  (§) 

A(7acid(l) 

AGacid  (1) 

acid 

(kJ/mol) 

SMD-GIL  (kJ/mol) 

SMD-H2O  (kJ/mol) 

HNCS 

1332.5 

HNCS 

607.0 

HNCS 

549.5 

HNO3 

1306.7 

HNCNCN 

578.0 

HNCNCN 

521.8 

HSCN 

1274.2 

HNO3 

566.6 

Hvdca  04 

515.2 

HNCNCN 

1272.6 

Hvdca  01 

565.2 

Hvdca  01 

514.6 

NCNHCN 

1234.6 

Hvdca  04 

564.9 

NCNHCN 

495.9 

Hvdca  01 

1212.1 

NCNHCN 

550.3 

HNO3 

494.9 

Hvdca  04 

1203.3 

HSCN 

548.8 

HSCN 

492.6 

HTCM  (central) 

1197.1 

HTCM  (central) 

537.9 

HTCM  (central) 

490.5 

HTCM  (terminal) 

1188.7 

HTCM  (terminal) 

533.9 

HTCM  (terminal) 

483.9 

Hvdca  02 

1152.7 

Hvdca  03 

518.1 

Hvdca  03 

471.4 

Hvdca  03 

1146.1 

Hvdca  02 

517.6 

Hvdca  02 

468.6 

Hvdca  01  = 

terminal  NH,  Hvc 

ca  02  =  -CCNiH,  Hvdca  03  =  -CCN2H,  Hvdca  04 

=  central  NH 
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It  is  well  known  that  carbenes  are  highly  reactive  species, and  it  is  likely  that 

these  carbenes  formed  upon  heating  of  the  ILs  are  also  involved  in  the  formation  of  the  species 

detected  at  m/z  C^A'  -  27  amu  seen  in  the  dca'  and  SCN'  systems.  Additional  evidence  for  these 

species  was  seen  in  the  TGA-MS  of  EMIM^SCN'  (SI)  and  in  the  GC/MS  trace  of  BMIM^dca' 

(Figure  4d),  indicating  that  these  products  are  thermodynamically  stable.  By  determining  the 

experimental  VUV  ionization  potential  (IP)  of  the  m/z  27  photoion  to  be  =  13.4  ±  0.2  eV,  this 

mass  can  be  identified  as  HCN  (I.P.  =  13.4  eV^*).  Detailed  ab  initio  analyses  of  the  possible 

mechanisms  of  formation  of  m/z  142  +  HCN  form  EMIM^SCN'  are  shown  in  Figure  7.  The 

proposed  mechanism  first  involves  the  proton  transfer  to  the  SCN'  anion  to  form  HNCS,  which  is 

more  basic  than  HSCN  (Table  4).  This  is  followed  by  the  rapid  nucleophilic  attack  of  the 

carbene  at  the  carbon  in  HNCS.  The  S  can  then  migrate  to  C2  on  the  imidazole  ring,  followed 

by  HNC  elimination  from  C2,  which  rapidly  converts  to  HCN.  This  process  effectively  replaces 

the  C2  proton  in  the  ion  pair  with  a  sulfur  atom,  resulting  in  a  neutral  thioimidazole,  which  could 

explain  the  lack  of  m/z  1 10  for  EMIM^SCN'  in  the  VUV-TOFMS  spectra  if  the  reaction  is 

sufficiently  fast  once  the  carbene  is  formed.  This  scheme  has  been  generalized  to  CN-containing 

anions  in  Figure  8.  The  more  basic  the  anion  is,  and  to  a  lesser  extent,  the  more  electrophilic  the 

protonated  anion  is,  the  more  likely  this  reaction  is  going  to  occur.  These  species  are  detected  at 

temperatures  slightly  higher  than  the  detection  temperature  of  the  carbenes,  which  is  consistent 

with  the  small  activation  barrier  involved  in  the  nucleophilic  attack  of  the  carbene  on  the 

protonated  anion  (step  2  to  3  in  Figure  7).  This  type  of  mechanism  may  also  explain  the  addition 

of  an  oxygen  atom  to  the  EMIM  cation  in  the  thermal  decomposition  of  EMIM^Os’.^^ 

Although  sulfur  addition  to  C2  on  the  dialkylimidazolium  cation  through  an  n-heterocyclic 

carbene  (NHC)  intermediate  has  been  demonstrated  previously,"^^  to  our  knowledge,  this  is  the 
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first  report  of  the  NCN-ylidene  being  formed.  This  route  eould  possibly  be  useful  for  a  one-pot 
synthesis  for  these  NCN  imidazolium-ylidene  speeies  without  the  need  for  flammable  solvents. 


58.4 


Figure  7.  M06/6-31+G(d,p)  free  energy  reaction  profile  of  the  thermal  decomposition  of 
EMIM^SCN'  leading  to  the  formation  of  m/z  142  and  HCN.  Values  are  in  kJ/mol 


HN  N 

\ — - 

Figure  8.  Generalized  scheme  for  the  addition-elimination  reaction  through  a  carbene 
intermediate. 
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At  even  higher  temperatures,  the  third  major  decomposition  pathway  indicated  in  the 
VUV-PI-TOFMS  experiments  is  the  Sn2  alkyl  abstraction  by  the  anion,  as  has  been  proposed 
previously  in  the  1 -ethyl- 3 -methylimidazolium  bromide*^  and  1 -ethyl-3 -methylimidazolium 
nitrate^^  ionic  liquids.  Due  to  steric  constraints,  the  abstraction  of  the  methyl  group  on  the 
imidazolium  cation  is  expected  to  be  more  facile  than  the  bulkier  ethyl  or  butyl  groups,  and  this 
is  evidenced  by  the  lower  activation  energies  for  the  formation  of  ethylimidazole  or 
butylimidazole  versus  methylimidazole  (Table  2  and  supported  by  DFT  calculations).  Further 
evidence  that  this  Sn2  reaction  occurs  is  the  detection  of  CH3NCS  and  CH3CH2NCS  in  the 
thermal  decomposition  of  EMIM^SCN':  The  experimental  I.P.s  of  mass  73  and  87  are  8.2  ±  0.2 
eV  and  9.1  ±  0.2  eV,  respectively,  while  the  M06  I.P.  values  are  9.1  ±  0.2  eV  and  9.0  ±  0.2  eV 
for  R-NCS  (R=CH3,  CH3CH2-)  and  9.9  ±  0.2  eV  and  9.7  ±  0.2  eV  for  R-SCN  (R=CH3,  CH3CH2- 
),  giving  a  better  match  with  the  RNCS  species.  The  temperature-jump  experiments  described 
below  support  this  conclusion  as  well. 

The  appearance  of  the  m/z  137  peak  in  the  BMIM^  ILs  indicates  that  the  BMIM^  cation 
loses  2  hydrogens  quite  easily.  It  is  possible  that  the  mechanism  for  this  involves  the  formation 
of  the  BMIM:  carbene  species  which  could  rapidly  react  to  lose  a  second  hydrogen.  In  Figure  5, 
the  progression  from  the  intact  BMIM"^  cation  (m/z  139,  formed  from  dissociative 
photoionization  of  the  vaporized  ion  pair)  through  the  carbene  (m/z  138)  to  m/z  137  as  the 
temperature  is  raised  from  438  to  483  K  can  be  seen.  At  438  and  453  K,  the  predominant  peak  is 
at  m/z  139,  with  smaller  peaks  at  m/z  138  and  137.  At  468  K,  the  m/z  139  and  137  peaks  are 
nearly  equal  in  intensity,  but  the  ratio  of  the  m/z  139/138  peaks  does  not  change  significantly. 
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This  seems  to  indicate  that  if  the  m/z  137  species  is  formed  via  the  BMIM:  carbene  (m/z  138), 
the  carbene  reacts  quickly  to  form  m/z  137  and  is  not  detected.  At  483  K,  the  m/z  139  peak 
nearly  disappears  and  m/z  137  dominates,  indicating  that  the  rate  of  conversion  of  m/z  139  to 
137  has  surpassed  evaporation,  and  the  intensities  of  m/z  138  and  139  are  consistent  with  C 
isotopes  resulting  from  the  m/z  137  species.  An  interesting  possibility  is  that  the  length  of  the 
butyl  group  allows  the  end  of  the  butyl  group  to  access  the  C2  carbene  structure  via  rotations 
around  the  N-C  and  C-C  bonds  in  the  butyl  group  (see  Scheme  1),  whereas  in  the  EMIM:,  the 
ethyl  group  is  too  short  to  react  with  the  C2  carbene  and  does  not  lose  a  second  hydrogen: 


( 


Scheme  1. 

Bicyclic  ionic  liquids  have  been  prepared  from  the  BMIM  cation-based  ionic  liquids, which 
means  that  cyclization  of  the  butyl  group  should  be  kinetically  accessible  to  the  BMIM  carbene 
cation.  It  could  be  possible  that  m/z  137  is  the  result  of  conversion  of  the  BMIM"^  cation  to  the 
bicyclic  cation  species  (B^,  m/z  137)  in  the  condensed  phase  through  the  carbene  intermediate, 
which  can  then  evaporate  as  a  B'^A'  ion  pair  that  undergoes  dissociative  photoionization  to 
produce  m/z  137  and  A*.  If  the  resulting  B^A’  has  a  lower  enthalpy  of  vaporization  than 
BMIM'^A’  then  the  m/z  137  peak  would  dominate  the  mass  spectrum  versus  m/z  139.  Another 
possibility  for  the  formation  of  m/z  137  is  as  an  ionization  fragment  of  m/z  178  (see  GC/MS 
spectrum,  Figure  4d).  However,  additional  research  is  needed  to  understand  and  explain  this 
phenomenon. 

Temperature-jump  Fourier  Transform  Infrared  Spectroscopy 
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Additional  evidence  for  the  Sn2  alkyl  abstractions  is  seen  in  the  T-jump  experiment 

involving  BMIM^SCN',  which  in  Figure  6a  clearly  shows  the  formation  of  CH3NCS.  For 

BMIlVf^SCN',  the  intense  peak  at  2080  cm  '  and  smaller  peak  at  ~2200  cm  '  in  Figure  6a 

suggests  the  formation  of  SCN-CH3,  whereas  NCS-CH3  would  have  only  a  single  peak  at  2180 

cm'*/"'  which  is  consistent  with  the  VUV-PI-TOFMS  results  for  BMIM^SCN'.  No  HNCS  or 

HSCN  was  detected.  These  experiments  also  indicate  the  demethylation  of  the  imidazolium 

cation  to  form  butylimidazole  and  a  methylated  anionic  species  for  BMIM^TCM'  and 

BMIlS/f^dca'  ILs.  The  methylated  dicyanamide  (with  two  -CN  groups)  and  methylated 

tricyanomethanide  (three  -CN  groups)  products  likely  polymerize  when  heated  slowly,  which 

would  explain  their  absence  in  the  VUV-TOFMS  experiments,  while  methylated  thiocyanate 

(with  only  one  -CN  group)  cannot  polymerize.  For  BMIM^TCM',  Figures  6b  and  6c 

demonstrate  a  good  match  from  2100  to  2300  cm'  with  data  from  Finnerty  et.  al.  (used  with 

permission),  indicating  methylation  of  the  terminal  nitrogen  on  the  tricyanomethanide  species, 

CH3-NCC(CN)2  rather  than  at  the  central  carbon,  CH3-C(CN)3.  This  is  consistent  with  the  M06 

IR  simulation  in  Figure  6e  indicating  that  there  should  be  no  peaks  in  this  region  for  the  centrally 

methylated  species.  For  BMIM^dca'  (Figure  6d),  there  are  two  peaks  which  arise  from  the 

methylated  product  at  ~1600  cm'*  and  ~2200  cm'*.  From  the  M06  IR  simulation  (Figure  6f), 

both  the  terminally  and  centrally  methylated  dicyanamide  species  show  peaks  in  the  2200  cm'* 

range,  so  the  indentity  of  this  product  isn’t  clear.  However,  the  simulated  spectrum  of  the 

terminally  methylated  species  shows  a  small  peak  at  ~1600  cm  '  whereas  the  centrally 

methylated  dicyanamide  simulation  does  not  show  any  peaks  in  this  region.  Also,  previous  DFT 

calculations  indicate  that  the  barrier  to  form  the  terminally  methylated  dicyanamide  should  be 

lower  than  for  the  centrally  methylated  product.^^ 
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Density  Functional  Theory  Modeling 

In  order  to  verify  the  thermal  deeomposition  meehanisms  proposed  based  on  the 
experimental  evidenee,  hybrid  density  functional  theory  was  applied  at  the  M06/6-31+G(d,p) 
level  of  theory  to  these  IL  systems.  M05  has  been  shown  to  provide  a  good  level  of  theory  for 
calculating  energetics  of  ionic  liquid  systems  (mean  average  deviation,  MAD,  of  9.5  kJ/mol),*^ 
and  M06  is  presumed  to  be  as  good  if  not  better  than  M05  for  our  purposes.  Stationary  states  of 
reactants,  intermediates,  transition  states  and  products,  as  well  as  adiabatic  photoionization 
potentials  were  calculated.  A  representative  free  energy  profile  for  EMIM'^SCN'  is  shown  in 
Figure  7. 

Adiabatic  ionization  potentials  (I.P.)  calculated  for  the  proposed  species  match  well  with 
the  experimentally  determined  photoionization  energies  (Table  3).  For  FMIM:  (m/z  1 10),  the 
experimentally  determined  value  is  7.6  ±  0.2  eV  and  the  M06  value  is  7.7  ±  0.2  eV.  Similarly, 
for  the  addition-elimination  reactions  (S:  m/z  142  and  170,  -NCN:  150  and  178  for  FMIM"^  and 
BMIM^,  respectively,  Figure  8),  all  of  the  experimental  I.P.s  match  the  M06  calculated  I.P.s 
within  0.2  eV.  In  some  cases,  these  species  are  detected  at  energies  below  the  threshold  to  detect 
from  dissociative  ionization  of  the  ion  pair,  which  is  an  indication  that  these  species  are 
formed  in  the  condensed  phase  and  are  not  a  result  of  dissociative  photoionization. 

Polarizable  Continuum  Model:  SMD-GIL  Method 

In  order  to  accurately  represent  the  energetics  involved  in  the  liquid  phase  reactions  of 

ionic  liquids,  a  comparison  was  made  between  the  gas  phase  AGacid(g)  calculations  (Table  4)  and 

the  M06  energetics  in  Figure  7a  and  the  corresponding  energetics  using  the  recent  adaptation  of 

the  continuum  solvation  model  (SMD)  referred  to  as  the  generic  ionic  liquid  (SMD-GIL) 

model.  This  model  parameterizes  the  PCM  to  accurately  determine  free  energies  of  solvation 
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and  benchmarks  to  known  experimental  values,  resulting  in  mean  unsigned  error  (MUE)  in  free 
energy  values  of  less  than  1.85  kcal/mol.  For  EMIM^dca',  the  SMD-GIL  solvent  model  has  a 
MUE  of  0.65  kcal/mol.  For  determining  the  AGacid(l)  in  the  condensed  phase,  the  estimated  error 
in  these  calculations  is  ±V3(0.65)  kcal/mol,  or  ±1.1  kcal/mol  (±4.7  kJ/mol).  The  results  for 
AGacid  are  seen  in  Table  4. 

The  results  of  the  SMD-GIL  calculations  for  EMIM^SCN’  are  seen  in  Figure  9.  It  is 
interesting  to  note  that,  relative  to  the  gas  phase,  the  free  energy  of  solvation  using  the  PCM 
model  stabilizes  the  ion  pair  much  more  than  the  other  neutral/molecular  species  and  transition 
states.  Intuitively,  this  makes  sense,  as  the  ionic  solvent  will  more  readily  solvate  ions  versus 
molecular  species.  The  effect  of  lowering  the  free  energy  of  the  IE  ion  pair  reactant  more  than 
the  free  energies  of  the  intermediates,  transition  states  and  products  increases  the  initial  free 
energy  of  activation  in  the  condensed  phase  relative  to  the  gas  phase.  Interestingly,  the  gas  phase 
prediction  of  the  activation  barrier  for  the  Sn2  reactions  for  EMIM^SCN'  incorrectly  predicts  the 
methyl  abstraction  to  have  a  higher  barrier  than  the  ethyl  abstraction,  but  the  PCM  model 
correctly  predicts  the  methyl  abstraction  to  be  the  more  facile  route. 
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^  HNC  - ►  HCIM 


EMIM^SCN-  ^  EIM  +  MeNCS,  AGt  =  148.9  (160.0) 
EMIM^SCN-  ^  MIM  +  EtNCS  ,  AQt  =  143.0  (166.5) 


Figure  9.  Free  energy  reaction  profile  comparison  of  gas  phase  (black)  and  SMD-GIL  (red) 
using  M06/6-3 1+G(d,p)  for  the  formation  of  m/z  142  from  EMIM^SCN'.  Values  are  in  kJ/mol. 


Nucleophilicity  of  Anions 

While  nucleophilicity  scales^^  are  well  defined  for  many  species  in  common  solvents  in 
which  synthetic  reactions  are  carried  out  (i.e.,  in  H2O,  hexane,  etc.),  the  relative  nucleophilicity 
of  the  nitrile  containing  anions  towards  alkyl  abstraction  from  the  dialkylimidazolium  cations  is 
not  readily  apparent.  In  order  to  evaluate  the  nucleophilicity  of  these  CN  containing  anions, 
transition  states  of  the  Sn2  type  methyl  abstractions-the  major  alkyl  abstraction  reaction  from  the 
dialkylimidazolium  cations  (Scheme  2)-  were  located  in  both  the  gas  phase  and  in  the  condensed 
phase  using  SMD-GIL  to  determine  their  free  energies  of  activation  (Table  5).  The  pathways 
with  the  lowest  activation  free  energies  would  be  expected  to  have  the  fastest  alkyl  abstraction 

rates  and  therefore  be  the  most  nucleophilic.  The  free  energies  of  activation  for  Sn2  methyl 
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abstraction  for  SCN',  dca',  TCM'  and  Br  (for  comparison  to  our  previous  study)  to  form 
CH3NCS  CH3N=C=NCN,  CH3N=C=C(CN)2  and  CH3Br,  respectively,  are  seen  in  Table  5. 


X  =  Br,  SCN,  dca,  TCM 


Scheme  2. 


Table  5.  Free  energies  of  activation  for  Sn2  methyl  abstraction  from  EMIM"^  by  the  anion 
calculated  at  the  M06/6-31+G(d,p)  level  of  theory  and  using  the  generic  ionic  liquid  (GIL)  model 
to  simulate  the  condensed  phase  methyl  abstraction. 


IL 

AG*  (gas:  M06)  (kJ/mol) 

AG*  (liquid:  GIL)  (kJ/mol) 

AAG  (gas^liquid)  (kJ/mol) 

EMIM^Br' 

138.6 

134.6 

-4.0 

EMIM^dca' 

163.6 

158.2 

-5.4 

EMIM^SCN' 

148.9 

160.0 

11.1 

EMIM^TCM' 

175.8 

181.5 

5.7 

From  Table  5,  the  trend  of  relative  nucleophilicities  of  the  anions  is,  in  decreasing  nucleophilic 
strength:  Br  >  SCN'  >  dca'  >  TCM'.  It  is  interesting  to  note  that  difference  in  the  free  energy  of 
activation  does  not  change  significantly  (from  -4.0  to  11.1  kJ/mol)  upon  going  from  the  M06  gas 
phase  calculation  to  the  GIL  condensed  phase  calculation,  which  contrasts  with  the  large  shift  in 
the  free  energy  of  activation  for  the  initial  proton  transfer  in  the  SCN'  system  (Figure  9).  Also, 
the  trend  for  anion  basicity  and  relative  nucleophilicity  are  the  same,  but  it  is  not  apparent  if  this 
is  coincidental  or  if  this  is  a  general  trend.  While  the  GIL  model  is  an  improvement  over  the  gas 
phase  approach  to  calculating  the  relative  anion  nucleophilicities  proposed  here,  the  gas  phase 
calculations  predict  the  same  trend  as  the  GIL  method,  and  could  be  considered  a  good  first 
approximation  to  determining  relative  anion  nucleophilicities. 

CONCLUSIONS 
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The  thermal  deeomposition  of  alkylimidazolium  ionic  liquids  with  CN  containing  anions 
has  been  investigated  by  multiple,  complementary  experimental  techniques  and  DFT 
calculations.  Two  common  thermal  decomposition  mechanisms  have  been  identified-carbene 
formation  and  alkyl  abstraction-  that  are  related  to  anion  basicity  and  nucleophilicity, 
respectively.  A  third  mechanism,  that  likely  proceeds  through  the  carbene,  allows  for  addition  of 
S  or  -NCN  groups  to  C2  on  the  imidazole  ring.  M06  calculations  in  the  gas  phase  and  the  SMD- 
GIL  variant  of  the  polarizable  continuum  model  for  the  condensed  phase  have  demonstrated  the 
ability  to  predict  trends  in  anion  basicity  and  anion  nucleophilicity  which  should  prove  useful  in 
predicting  thermal  stability  trends  in  dialkylimidazolium  ionic  liquids  and  perhaps  could  be  used 
to  predict  thermal  stabilities  of  ionic  liquids  in  general. 
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Supplemental  Information 


BMIM^dca 


temperature  (K) 


Figure  SI.  TGA  measurements  for  BMIM^dea':  a)  non-isothermal  TGA,  b)  A//(overall)  from 
(a),  e)  A//vap  from  isothermal  TGA. 
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Figure  S2.  a)  TGA-MS  of  m/z  97  and  82,  and  b)  Af/^  of  product  formation,  and  c)  complete 
TGA-MS  at  493,  518  and  543  K. 
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Figure  S3.  VUV-PI-TOFMS  of  BMIM^dca  a)  at  438  K,  b)  an  expanded  view  of  (a),  c)  at  483 
K,  d)  an  expanded  view  of  (c),  and  e)  temperature  profdes  of  photoion  peak  intensities. 


Figure  S4.  T-jump  experiment  of  BMIM^dca'. 
EMMIM^dca’ 


t«mperature  (K) 


Figure  S5.  TGA  measurements  for  EMIM^dca':  a)  non-isothermal  TGA,  b)  A//(overall)  from 
(a),  c)  AHyap  from  isothermal  TGA. 
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Figure  S6.  a)  TGA-MS  of  m/z  96  and  110,  and  b)  of  product  formation  in  (a). 


Figure  S7.  VUV-Pl-TOFMS  of  EMMlM+dca  a)  at  423  Kand  10.0  eV,  b)  at  498  K  and  10.0  eV, 
and  c)  temperature  profiles  of  photoion  peak  intensities. 
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Figure  S8.  TGA  measurements  for  EMEVt^SCN’:  a)  non-isothermal  TGA,  b)  A//(overall)  from 
(a),  e)  A//vap  from  isothermal  TGA. 
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Figure  S9.  a)  TGA-MS  of  m/z  73  and  87,  b)  A/7^  of  product  formation  in  (a),  and  c)  complete 


TGA-MS  at  473,  498  and  523  K. 


Figure  SIO.  VUV-PI-TOFMS  of  EMIM+SCN'  at  467  K  and  10.0  eV. 
BMIM+SCN' 


Figure  Sll.  TGA  measurements  for  BMEVt^SCN':  a)  non-isothermal  TGA,  b)  A//(overall)  from 
(a),  c)  AHvap  from  isothermal  TGA  and  d)  TGA-MS. 
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Figure  S12.  VUV-PI-TOFMS  of  BMIM+SCN'  at  499  K  and  10.0  eV. 
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BMIM+TCM' 
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Figure  S13.  TGA  measurements  for  BMIM^TCM’:  a)  non-isothermal  TGA,  b)  A//(overall) 
from  (a),  c)  isothermal  TGA,  and  d)  Affyap  from  isothermal  TGA.  (reproduced  from  J  Phys  Chem 
Lett,  used  with  permission). 
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Figure  S14.  a)  TGA-MS  of  m/z  78  and  82,  b)  AH^  of  product  formation  in  (a),  and  c)  complete 
TGA-MS  at  498,  523  and  548  K. 


Figure  S15.  VUV-PI-TOFMS  of  BMIM^TCNT  at  433  K  and  10.0  eV. 
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Figure  S16.  TGA  measurements  for  BMIM\dca‘:  a)  non-isothermal  TGA,  b)  A//(overall)  from 
(a). 
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Figure  S18.  VUV-PI-TOFMS  of  BMIMVdca  at  523  K  and  10.0  eV. 
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